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Abstract 



We study the viscosity corrections to the growth rate of nucleating bubbles 
in a slightly supercooled first order phase transition in 1 + 1 and 3 + 1 di- 
mensional scalar field theory. We propose a microscopic approach that leads 
to the non-equilibrium equation of motion of the coordinate that describes 
small departures from the critical bubble and allows to extract the growth 
rate consistently in a weak coupling expansion and in the thin wall limit. 
Viscosity effects arise from the interaction of this coordinate with the stable 
quantum and thermal fluctuations around a critical bubble. In the case of 
1 + 1 dimensions we provide an estimate for the growth rate that depends on 
the details of the free energy functional. In 3 + 1 dimensions we recognize 
robust features that are a direct consequence of the thin wall approximation 
that trascend a particular model. These are long-wavelength hydrodynamic 
fluctuations that describe surface waves. We identify these low energy modes 
with quasi-Goldstone modes which are related to surface waves on interfaces 
in phase ordered Ising-like systems. In the thin wall limit the coupling of 
this coordinate to these hydrodynamic modes results in the largest contri- 
bution to the viscosity corrections to the growth rate. For a </> 4 scalar field 
theory at temperature T < T c the growth rate to lowest order in the quartic 



with R c ; £ the critical ra- 



self-coupling A is O = ^ 1 - 0.003AT£ 

dius and the width of the bubble wall respectively. We obtain the effective 
non-Markovian Langevin equation for the radial coordinate and derive the 
generalized fluctuation dissipation relation, the noise is correlated on time 
scales il™ 1 as a result of the coupling to the slow hydrodynamic modes. We 
discuss the applicability of our results to describe the growth rate of hadron 
bubbles in a quark-hadron first order transition. 
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I. INTRODUCTION 



The dynamics of first order phase transitions is a fundamental ingredient in particle 
physics and in condensed matter. First order phase transitions occur via the nucleation of 
bubbles of the true vacuum phase in the metastable or false vacuum phase. At large tempera- 
ture it is mediated by thermal or overbarrier activation and at low temperatures via quantum 
nucleation. First order phase transitions are conjectured to occur in QCD and in the Elec- 
troweak theory. In QCD a first order phase transition could describe the hadronization of the 
quark-gluon plasma, possibly produced in the early Universe at about 10~ 5 seconds after the 
Big Bang or in relativistic heavy ion collisions |]T]-§| • In the Electroweak theory a first order 
phase transition is argued to provide the non-equilibrium setting for baryogenesis |§-f|. In 
early universe cosmology, first order phase transitions had been proposed as mechanisms to 
generate the inflationary stage |7]-|9| • In condensed matter physics thermal activation results 
in the nucleation of bubbles of the lowest free energy phase in binary fluids pH|jn|] and also 
of the decay of metastable dimerized states in quasi-one dimensional charge density wave 
systems and non-degenerate organic conductors fll2,13|. 

The most comprehensive microscopic theory of nucleation via thermal activation was 
presented by Langer [O] and was later extended to quantum field theory to account both 



for thermal activation as well as for quantum nucleation ||To]-pifj. An approach to describe 
nucleation in a non-steady state situation in real time has been advocated in . In the limit 
in which nucleation is dominated by thermal activation and overbarrier transitions, these 
bubbles are produced via large thermal fluctuations. These bubbles will grow whenever their 
radius is larger than a critical value and collapse if it is smaller. Supercritical bubbles will 
grow to convert the metastable phase into the stable phase or until they percolate achieving 
the total conversion of the metastable phase. For slightly supercritical bubbles an important 
dynamical quantity is the growth rate of a bubble 



dt 



In 



R(t) 



R(t) > R c 



with R(t) the (time dependent) radius of a slightly supercritical nucleating bubble and R c 
is the critical radius [l4]||,|llj. Langer's theory provides the nucleation rate per unit volume 
per unit time given by 



[1.2) 



where Ft, is the free energy of a critical bubble and T> is proportional to the ratio of the 
determinants of the fluctuation operators around the bubble configuration to that around 



the homogeneous metastable state ||T^,|3|JT0|| . The regime of applicability of homogeneous 
nucleation theory is for Fb » T, for Fb ~ T small amplitude thermal fluctuations can 
trigger the phase transition and nucleation and spinodal decomposition can no longer be 
distinguished. The ratio of determinants in T> has been computed analytically and numer- 
ically in several important cases [f2C| -|2"2|l . Csernai and Kapusta have studied the growth 
rate of hadronic bubbles in a quark gas by extending Langer's theory to the relativistic case. 
These authors studied a coarse-grained field theory in terms of a relativistic hydrodynamic 
description with viscous terms in the energy momentum tensor and the baryon current. 
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Their conclusion was that the growth rate is determined by the coefficients of the shear 
and bulk viscosities and that hadronic bubbles do not grow when these coefficients vanish. 
Another approach presented by Ruggeri and Friedman [^3[] also based on baryon-free rel- 
ativistic hydrodynamics but with a different treatment of the heat conduction and energy 
flow reached a different conclusion, th at the growth rate is non-zero even for vanishing bulk 
and shear viscosities and that viscosity effects are subleading for small viscosities. If the 
hadronization phase transition is of first order there are potentially important experimental 
signatures associated with the homogeneous nucleation of the hadronic phase some of which 
had been studied in reference 



Homogeneous nucleation of the quark-gluon plasma has been recently studied with a bag 
model of the equation of state for the quark and hadron phases and the different predictions 
above have been compared |25j. The formulation and results of Csernai and Kapusta had 
recently been used to study first order quark-hadron phase transition in the Early Universe 
for a first order hadronization transition and more recently homogeneous nucleation 
has been tested numerically in 2 + I dimensional systems with qualitative agreement to the 
standard result ]p7[|. 

An alternative phenomenological description of the growth rate based on dissipative 
hydrodynamics combined with the finite temperature effective potential has been presented 
in [28]. In this work analytical and numerical studies reveal a dependence of the growth rate 
on the phenomenological dissipative coefficient. 

Khlebnikov |29| and Arnold |30] have studied dissipative effects that slow down the 



growth rate of a supercritical bubble by coupling the order parameter that describes nu- 
cleation to other fields (with a trilinear coupling) and applie d the fluctuation dissipation 
relation to one loop order in the coupling to the other field. 

Finally we must mention a numerical approach to the description of nucleation, based on 
a phenomenological Markovian Langevin dynamics in terms of the finite temperature effec- 
tive potential (typically computed to one loop order) with a (local) friction coefficient and 



a Gaussian white noise term related by the fluctuation-dissipation theorem [31,32]. The use 



of the finite temperature effective potential in the description of the inhomogeneous bubble 
configuration that seeds the nucleation of the phase of lowest free energy is an important 
ingredient, in particular in coarse-grained phenomenological descriptions Q. Finite temper- 
ature effects are included in the coarse grained free energy which describes long wavelength 
physics by integrating o ut short wavelength fluctuations which are in thermodynamic equi- 
librium [P3| , |3"3{] . These finite temperature corrections modify quantitatively and sometimes 
substantially the equilibrium free energy functional, for example the position of the min- 
ima, masses and couplings. These are important parameters that determine the profile of 
the bubble configuration, since the asymptotic behavior of the bubble is determined by the 
position of the minima of the equilibrium free energy. 

The importance of the growth rate for describing the dynamics of nucleation in the 
hadronization and the Electroweak phase transitions as well as the practical importance 
of nucleation in quasi-one dimensional organic conductors justifies a continued effort to 
understand from a microscopic perspective the influence of dissipation, viscosity and friction 
upon the growth rate of supercritical nucleating bubbles. 

Focus and Strategy: Our goal in this article is precisely to study dissipative effects 
upon the growth rate Q from a more microscopic point of view in model quantum field the- 
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ories that describe the main features of nucleation. The simplest model to study nucleation 
and a first order phase transition is a 4 field theory with a small explicit symmetry breaking 
term that breaks the degeneracy between ground states and thus leads to the existence of a 
metastable state. 

Although arguably this model could hardly describe the features of the quark-hadron or 
electroweak phase transitions, our hope is to extract robust phenomena that will be generic 
to the physics of nucleation and that could enlighten the effect of viscosity on the growth 
rate. 

The study begins by identifying the inhomogenous bubble configuration which is a so- 
lution of the static equations of motion. We include the finite temperature effects by con- 
sidering the equation of motion in terms of the finite temperature effective potential, this is 
achieved consistently by the addition of finite temperature counterterms to the Lagrangian. 
The quadratic fluctuations around this bubble solution feature an unstable mode that de- 
scribes small perturbations from the critical bubble. The instability is a manifestation of 
the growth (or collapse) of supercritical (or subcritical) bubbles and the unstable eigenvalue 
is directly related to the growth rate of slightly supercritical bubbles. Besides this unstable 
mode there are modes of zero frequency corresponding to translations of this configuration 
and modes of positive frequency that correspond to stable fluctuations. 

The main strategy is to consider the dynamics of the coordinate that determines the 
unstable direction in functional space, and treating it as a collective coordinate. We then 
focus on the interaction of this coordinate with those representing the stable fluctuations by 
expanding the original Lagrangian in terms of these coordinates in function space and rec- 
ognize the interaction terms between these coordinates. The translational mode is anchored 
and ignored since we are not interested in studying the dynamics of the translational degree 
of freedom but that of the growth of a supercritical bubble. We assume that the stable 
degrees of freedom are in thermal equilibrium and that they serve as a "bath" with which 
the unstable coordinate interacts. Integrating out the degrees of freedom associated with 
the stable fluctuations we obtain an effective description of the dynamics for the unstable 
coordinate which includes the "viscosity and friction" effects associated with the transfer of 
energy to the stable modes. This description allows us to obtain the equations of motion for 
the expectation value of the unstable coordinate wherefrom we can obtain the corrections 
to the growth rate from the interaction with the stable degrees of freedom. Furthermore, 
integrating out the stable degrees of freedom in the non-equilibrium functional allows us to 
extract the effective Langevin equation for the unstable coordinate and to obtain consis- 
tently the noise terms. We emphasize that this approach is fundamentally different from the 
previous studies described above. First of all, we are not computing the ratio of determinants 
that enter in the rate expression. This ratio of determinants only involves the quadratic fluc- 
tuations but not the interaction between the coordinates associated with the fluctuations. 



Our approach is different from that of references [[29|,[$(J], in these references the coupling 
of the order parameter associated with the nucleating field was to other independent fields. 
Contrary to this we consider the coupling of the unstable coordinates to the stable fluctu- 
ations involving the same field. Our approach starts from the microscopic theory and does 
not rely on a hydrodynamic description, however as it will be seen later, long-wavelength 
hydrodynamic fluctuations associated with surface waves provide the largest contribution to 
viscosity effects in three spatial dimensions. 



4 



In a very specific sense our program obtains the effective action for the degree of freedom 
that represents departures from a critical bubble by integrating out the degrees of freedom 
describing stable fluctuations, but of the same field. Our approach must necessarily rely 
on several different approximations: i) weak coupling to allow a perturbative expansion of 
the effective action and the real time equations of motion, ii) the thin wall approximation 
in which the critical radius is much larger than the width of the bubble wall. This ap- 
proximation is necessary to be able to provide quantitative answers. Finally as it will be 
justified later the important fluctuations can be described by the classical limit of the finite 
temperature distribution functions. 

We carry out this program to lowest order in the quartic coupling and to leading order 
in the thin wall approximation both in 1 + 1 and 3 + 1 dimensions. The motivation to study 
the lower dimensional case is provided by its potential application in experimentally relevant 
condensed matter systems |12||13l. 

Summary of results: 

Our main results can be summarized as a consistent perturbative correction to the growth 
rate from "dynamical viscosity" effects associated with the interaction of the radius of the 
bubble with the stable quantum and thermal fluctuations. In 1 + 1 dimensions the growth 
rate for bubbles is given by eqn. ( |3.15| ) with R c , £ the radius and width of the critical bubble 
respectively, A the quartic self-coupling and F[R c /£\ a slowly varying dimensionless function 
that depends on the details of the potential and spectrum of stable fluctuations, in particular 
scattering states. 

In 3+1 dimensions we find that in the thin wall limit the important low energy excitations 
are hydrodynamic fluctuations associated with surface waves. These are identified as quasi- 
Goldstone bosons and give the dominant contribution to viscosity effects on the growth rate. 
These low energy excitations are a robust feature of any model and not particular to the one 
considered in this article, however, the coupling of the unstable radial coordinate to these 
fluctuations depends on the model. To lowest order in the coupling the growth rate is given 
by eqn. flPOp . 

Furthermore, we obtain the effective Langevin equation for the unstable coordinate that 
reveals the correct microscopic correlations of the noise term showing that the coupling to 
the hydrodynamic modes determine that the noise is correlated on time scales Q^ 1 . 

The article is organized as follows: section II introduces the model, describes the strategy 
in detail and sets up the general form of the correction to the growth rate to lowest order (one 
loop) and in the classical limit. In section III the case of 1 + 1 dimensions is analyzed in detail 
within the model considered. In section IV we study the 3 + 1 dimensional case, and discuss 
the low energy fluctuations associated with surface waves, providing the argument that these 
are quasi- Goldstone hydrodynamic fluctuations. We argue that these low energy fluctuations 
are present independent of the model and will dominate the viscosity corrections to the 
growth rate. The conclusions are presented in section V, along with a critical discussion 
on the validity of our results to describe the growth rate of supercritical bubbles in coarse- 
grained descriptions of a first order quark-hadron phase transition. An appendix presents 
the effective Langevin equation for small departures of the critical radius and analyzes the 
generalized fluctuation-dissipation relation. 
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II. THE MODEL: GENERALITIES AND STRATEGY 



We consider a real scalar field, 4>(x), whose dynamics is determined by the following 
Lagrangian density 

£ = ^(^0)(^0)-y(0), (2.i) 

where V((f>) is a double well potential with a metastable minimum at 0_ and a stable 
minimum at + . It proves convenient to parametrize it in the following form: 

V(0)=A(0-0_) 2 0(0-0*). (2.2) 

This form of the potential is depicted in fig.([|) and is very similar to the coarse-grained 
free energy proposed in || identifying the local energy variable e in that reference with 
the scalar field 0, and also similar to the effective potential description used in numerical 



simulations for the Electroweak theory pl| , |32| . |34 1 . 
The stable minimum + is given by 

1 



+ 30, - ^402 + 90 2 - 40_0*J . (2.3) 

where the point 0* is parametrized by the mass m of small oscillations around the metastable 
minimum </>_ as 



where 



e = 



9 

m 



2A0 2 



;i - c) (2-4) 



and m 2 = V" (((>-) (2.5) 



and A is the quartic coupling. Although this parametrization does not look familiar it will 
prove advantageous later. In particular, (1 — e) gives a measure of the free energy difference 
between the true and the false vacua and is therefore related to the supercooling temperature. 
As e — > 1, the two minima become degenerate and the condition e ~ 1 defines the thin- wall 
limit (this will become clear later when we analyze the bubble profile), when the radius 
of the nucleating bubble is much larger than the width of its surface £ (or the correlation 
length). This corresponds to the case of small supercooling T < T c but T/T c pa 1 in the 
description of reference ||. In the thin wall limit, we find that + has the following simple 
form 

+ = -%( e -l) + 0(( e -l)2). (2.6) 
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A. The counterterms: 



We anticipate that there will be renormalization of the parameters in the potential, 
not only arising from ultraviolet divergences but more importantly from medium effects, 
i.e. finite temperature contributions to the mass, couplings and explicit symmetry breaking 
terms. In particular in three spatial dimensions we expect a correction to the mass and 
the linear symmetry breaking term both proportional to T 2 ,T, The origin of these 
finite temperature corrections are the usual tadpole diagrams, the correction to the linear 
symmetry breaking term is a consequence of the trilinear coupling. These finite temperature 
corrections are the usual ones obtained in an effective potential description and provide a 
finite renormalization to the potential. Our aim is to describe the dynamical aspects of the 
thermal fluctuations that are responsible for thermal activation and the nucleation of the 
true phase in the false vacuum phase. These are the solutions of the static field equations 
but in presence of the thermal fluctuations, i.e. with the potential renormalized by the finite 
temperature effects. That is to say, the potential that must enter in the field equations 
must be the finite temperature effective potential |33|j3l|j3^j27|j3^] . The finite temperature 



effective potential includes the finite temperature corrections for static and homogeneous 
field configurations. Including these finite temperature effects in the equations o f motion 
guarantees that the inhomogeneous bubble configurations will asymptotically tend to the 
homogeneous values of the field representing the correct extrema of the equilibrium free 
energy. Furthermore, these finite temperature renormalizations will determine the correct 
scales for the bubble size, its width and the unstable frequency associated with the growth 
of a supercritical bubble. 

We account for these finite temperature renormalizations by replacing the parameters in 
the potential by the finite temperature renormalized parameters and adding counterterms 
to the Lagrangian. The counterterms are then found consistently in a coupling or loop 
expansion by requesting that they cancel the contributions from the tadpole (or similar) 
terms completely, i.e. including the finite temperature parts. This prescription is at the 
heart of using the finite temperature effective potential to study the solutions that lead to 
the decay of the metastable state P J55||^ , p£l| . |2"7j , p2| . Thus the Lagrangian density becomes 



m = \(d^)(dy) - V{<f>,T) + 5CM (2.7) 



where the potential V(<f>,T) is of the same form as in (2.2) but in terms of the parameters 



that include the finite temperature (and ultraviolet) renormalizations, i.e. X(T) , m(T) , e(T) 
etc. and the counterterm Lagrangian density is of the form 

5C ct = + 5g(f) 3 + 5M 2 (f) 2 + H<f> (2.8) 

where space-time translational invariance dictates that the counterterms are constant. These 
counterterms will be consistently chosen in the loop expansion to cancel the contribution 
from local tadpoles. We will argue below that this procedure renormalizes consistently in 
perturbation theory the static properties of the bubble, in particular the width and radius of 
a critical bubble. In particular the free energy density studied in P, pT| , |32| , p4| , |33| is precisely 
of this form with parameters that depend on temperature. 
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At this point we may consider adding to the counterterms a wave-function renormal- 
ization which would involve corrections for inhomogeneous configurations. A wave-function 
renormalization will arise at one loop order and beyond because of the trilinear coupling. 
However, computing this correction in the usual loop expansion around a homogeneous 
background is not very relevant for the bubble configuration. Since our goal is to obtain the 
real time equations of motion for particular fluctuations around the bubble configuration, 
these renormalizations will arise automatically in the equations of motion obtained to one 
loop order. Thus the counterterms that we include in the Lagrangian density are only those 
for static renormalization thus accounting for the correct scales and asymptotic values of 
the classical solution and obtain the equivalent of the wave-function renormalization directly 
from the real time evolution of fluctuations around the bubble configuration. This procedure 
will become clear when we obtain explicitly the evolution equations in the next sections. 

We will focus our attention on the cases of one and three spatial dimensions. Although 
the case of one spatial dimension is of limited interest in particle physics, it is important 
in condensed matter and statistical physics of quasi-one dimensional systems |12| , |T3|| . Fur- 



thermore most aspects of the treatment are general and the one- dimensional case provides 
a somewhat simpler setting to introduce the main strategy as well as to explore the general 
features. The three dimensional case will be studied in detail subsequently. 



B. The bubbles: 

Before proceeding to the specific situations, we now consider general features of the 
bubble solutions and the quadratic fluctuations around them to focus on the precise strategy 
to follow. 

A critical bubble is a solution of the static field equations, which in terms of the (effective) 
potential V(4>, T) is given in one spatial dimension by 



d 2 (j) b dV{<j>,T) 



dx 2 



= (2.9) 

4>b 



with boundary conditions such that (fi b — > for \x\ — ► oo. Such a solution describes 
a "bubble" like configuration which approaches the false vacuum at asymptotically large 
distances and probes the true vacuum in a localized region in space of size 2R C with R c 
the critical "radius". Such a solution will be found in detail in section III. In three spatial 
dimensions, the critical bubble is a radially symmetric static field configuration solution of 
the following equation 



d 2 (j) b 2 d(f) b + dV(<j>,T) 



dr 2 r dr 



(2.10) 



with the boundary condition <p b (r) — ► 0_ as r — > oo. It corresponds to a field configuration 
that starts close to the true vacuum, <p + , and tends to the false vacuum, </>_, at asymptot- 
ically large radial distance. The solution of (|2.10 ) will be studied in detail in the thin wall 
approximation in section IV. 

In both cases the change from the true vacuum to the false one occurs around the ra- 
dius of the bubble, R c , over a distance £(T) ~ 1/m that defines the wall thickness of the 
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bubble and is related to the correlation length in the metastable phase. We will study the 
non-equilibrium dynamical viscosity of nucleating bubbles in the thin wall limit in which 
the radius of the bubble is much larger than the wall thickness; i.e. R 3> £. These field 
configurations are parametrized by two important coordinates: x which describes the posi- 
tion of the center of the bubble and the radius R, i.e. a bubble configuration is of the form 
4>(x — xq] R) and a critical bubble corresponds to R = R c determined by the solution to the 
equations of motion ( |2.9|j2.1U| ). The coordinate xq is associated with translational invariance 



and typically treated as a collective coordinate p5| , while the radius R determines the size 
of the bubble and will be treated also as a collective coordinate (see below). Since we are 
not concerned here with the dynamics of the translational degrees of freedom, but rather 
with that of the radius, we "clamp" the collective coordinate xq by fixing the center of the 
bubble at xo = 0. Integration over this collective coordinate results in the typical volume 
factor [114, 1^1 and is not relevant for our discussion. 



In d space dimensions, we can use the radius of the bubble R as a variational parameter 
and introduce the total energy of a bubble configuration with radius R by 



E var {R) = J d d x 



^(vu^R)) 2 + v(M^R)) 



(2.11) 



and in the limit R/£(T) » 1 (thin wall) has two main contributions: a volume contribution 
proportional to R d which is negative and arises from the region of the bubble that probes 
the true vacuum which has lower (free) energy, and a surface term which is positive and 
arises from the region of the bubble corresponding to the wall, includes gradient terms and is 
therefore proportional to the area of the surface of the wall, R d ~ l . In one spatial dimension 
the "surface" of the bubble corresponds to two points and the gradient energy saturates to 
a constant independent of the radius for large bubble radius. 

The general behavior of the total energy of a bubble configuration as a function of 
s = R/£ is depicted in fig. (^) for one spatial dimension and in fig. (^) for three spatial 
dimensions (see sections III and IV). The maximum of the energy function determines the 
value of the critical radius R c as discussed in detail in the next sections. Clearly a critical 
bubble is an unstable static solution of the equations of motion. For R < R c the surface 
energy term dominates and the bubble shrinks into the false vacuum phase, for R > R c 
the volume energy dominates and the bubble will grow as the gain in the volume energy is 
greater than the cost in elastic surface energy. 

A more clear description is obtained in the case of three spatial dimensions for a spheri- 
cally symmetric bubble. Let a be the surface tension (energy per unit area) of the spherical 
bubble configuration separating the metastable phase from the stable one. Then the total 
energy of the bubble is 

A-jr 

E var (R) = AFR 3 + AiraR 2 + ■ ■ ■ (2.12) 

3 

where AjF = \V(<f) + ,T) — V(0_,T)|, and the dots stand for corrections that are subleading 
in the thin wall approximation, see section IV (eqn. (|4.22|) for details. These (small) 
corrections will be neglected for the arguments presented below. The energy attains its 
maximum E c at the critical radius 

2a 

ft=^ (2.13) 
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and it is given by 



E c = E var (R c ) = —^- R l = 3 (AJF) 2 ' ( 2 - U ) 

Near the maximum of the energy function, it can be written as an expansion in terms of the 
(small) departures from the critical radius as 

E var (R) = E C - ho\R - R c ) 2 + ■■■ (2.15) 

where the frequency 

2 ^ E var {R) 



to 



dR 2 



?trt (2.16) 



is independent of the critical radius and AjF and as it will be seen in detail below, it is 
related to the growth rate of slightly supercritical bubbles. 



C. Fluctuations 



The study of the fluctuations around the classical bubble configuration begins by studying 
the spectrum of the fluctuation operator 



d 2 



dx 

-V 2 + V"((j) h {xiR c )) 



U n (x) = LO n U n (x) in ID 



U n (x) = uJA n {x) in 3D 



(2.17) 
(2.18) 



where the prime represents differentiation with respect to the field <fi. 

Taking a spatial derivative of the equation of motion satisfied by the bubble configuration 
it is straightforward to find that 



U (x,R c ) oc 



R c 

dx 



is a solution of the one dimensional eigenvalue problem (|2.17 ) with Uq = and that 

x d(f) b (r,R c ) 



Uo(x, R c ) oc 



\x\ 



dr 



(2.19) 



(2.20) 



are eigenf unctions of the three dimensional fluctuation operator with zero eigenvalue. These 
are the zero modes arising from translational invariance [14,l^,p3|. In one dimension the zero 
mode is an odd function of x since the bubble solution is even, therefore it has a node. Hence 



there must be another solution of the Schroedinger like operator ( 2.17|) that has no nodes 
and has a smaller and therefore negative eigenvalue. The form of the bubble solution in the 
thin wall approximation is that of a widely separated kink-antikink pair, as it will be shown 
in detail in the next section we find that the eigenfunction corresponding to the negative 
eigenvalue must be W_i oc d(j)b(x,R)/dR\n c . This can be understood simply from the fact 



10 



that the zero mode associated with the bubble solution is the antisymmetric combination 
of the zero modes associated with the individual kinks [35], hence the eigenfunction with 



lower eigenvalue must be the symmetric combination and the form of the bubble profile 
immediately leads to the conclusion that is the derivative with respect to the radius (see 
section III, eqn. ( |3.11| )). In the three dimensional case, the zero modes ( p.20| ) correspond 
to the angular momentum I = 1 spherical harmonics, therefore there must be an I = 
(spherically symmetric) solution of the Schroedinger-like eigenvalue problem ( |2.18| ) with 
a negative eigenvalue. In a later section we study in detail the bubble solution and the 
spectrum of fluctuations around it and conclude that the eigenfunction corresponding to the 
negative eigenvalue is W_i oc d4>b(r, R)/dR\n c in the thin wall approximation (see section IV, 
eqn. fll27D for U 000 ). 

Therefore the respective spectra of the fluctuation operators are 



U- 1 (£,R e ) = JN- 1 



dcj)(x, R) 



dR 



uj\ = -Q 2 ; Q 2 > ; U (x,R c 



R c 



U n>0 (x,R c ), ujI>0 



'No V<j>(x,Rc), ul = 
(2.21) 



with iV_i; iVo normalization factors (chosen real). The classical bubble solution corresponds 
to a saddle point in functional space, the mode U-\ determines an unstable direction. 

The field operator can now be expanded in the complete basis of eigenfunctions of 
( ETLlEOl) in either case and write in general 



(x,t) = (f) b (x -x ,R c ) + g_i(t)W_i(a? - x ) + q {t)U {x - x ) + Qn{t)Ua{x - x ), 



n>0 



(2.22) 



Whereas the bound state eigenfunctions are chosen real, the scattering states are chosen to 
satisfy the hermiticity condition U n (x) = W n (x); q n (t) = q*L n (t). 

The quanta associated with the coordinates q n> o for the stable modes will be referred 
generically as "mesons" , and the operators q n >o create and annihilate meson states. 

From ( |2.21| ) it is clear that the modes Uq ; U-\ correspond to shifts in the position of the 
bubble, xq and the radius R respectively. The collective coordinate treatment absorbs 
the zero mode Uq into a definition of a new quantum mechanical degree of freedom Xo(t) 
and expands the field operator in terms of the directions perpendicular to this zero mode, 



(f)(x , t) = 4> b (x - x (t), R c ) + q. 



-l(t)W_i(£ - X (t)) + J2 Qn(t)U n (x - X {t)) 

n>0 



(2.23) 



We are not interested in the dynamics of the translational collective coordinate, therefore 
we will "clamp" the position of the bubble and focus solely on the dynamics of the unstable 
coordinate g_i(t). Technically this is achieved by inserting a functional delta function for 
Xo(t) in the path integral, the corresponding Jacobian from the change of field variables to 
the collective coordinate [ 35 1 leads to the usual volume factor [|l4| , |r5f , hence in what follows 
we will set Xo(t) = 0. We note that for small amplitudes of the unstable mode we can write 
the expansion ( |2.23| ) as 



(x,t) w <j) b {x,R(t)) + ]T q n {t)U n {x) ■ R(t) = R C + q-i(t)JN. 



[2.24) 



n>0 
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that allows us to identify B(t) — R c = g_i(t)\/iVZi and the coordinate g_i is associated with 
small departures from the critical radius. 



D. The relation between SI 2 and lu 2 



At this stage we recognize that the fluctuation mode oc dd b (x, B c )/dB c is the 

(unstable) functional direction along which the bubble either grows or collapses, therefore the 
coordinate g_i (t) oc {R(t) — R c ) describes the small fluctuations around the critical bubble. 
An immediate question is: what is the relation between the frequency u 2 in (|2.15| , |2~T6| ) and 
Q 2 , the eigenvalue associated with U-\ ( |2.21| )? The answer to this important question is 
obtained by taking the second derivative of the total energy in eq. fl2.11|) with respect to 
the radius of the bubble R, which results in the following equation 



d 2 E var (R) 
dR 2 



d d r 



V 



dR 



V^ft.V 



'd 2 



dR 2 



+ 



' d 2 s b \ dv{s b 

dR 2 ) 86 



+ 



'ds b y d 2 v{s h 



dR 



(2.25) 



Integrating the first two terms in eq. (|2.25|) by parts using the boundary conditions for 
the bubble solution and eqs. ( |2 . 9| , |2 . 1 (J| , |2 . 1 7| j2 . 1 ^ ) , which are evaluated at the critical radius 



R c , we obtain the following relation 

d 2 E var (R) 



dR 2 



-LU 



Q 2 



(2.26) 



in terms of the normalization factor of the unstable mode (see ( [2.21| )) which will be found 
for the one and three dimensional cases separately below. We note that the relation (|2.24|) 
determines that in the absence of interactions with other degrees of freedom the growth rate 
of a supercritical bubble is given by 



d_ 

dt 



In 



B(t) 



B c 



1 



n 



(2.27) 



for small departures from the critical radius. 

Equation ( |2.26| ) is very useful since it relates the second derivative of the total classical 
variational energy E var (R) (the energy of the classical bubble as a function of its radius) to 
the eigenfrequency associated with the quantum unstable mode. Usually it is a difficult task 
to find the spectrum of the quantum fluctuation operator, and consequently the frequency 
of the unstable mode. The above relation can be used to find the value of the unstable 
frequency, or the normalization factor N_i if the frequencies are known. 



E. The strategy 

Our goal is to study the effects of quantum and thermal fluctuations upon the real time 
evolution of a bubble whose radius is slightly larger than critical. This is achieved by treating 
the coordinate q~i{t) oc (R(t) — B c ) as a collective coordinate, i.e. the radius of the bubble is 
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treated as a fully quantum mechanical variable interacting with the other degrees of freedom 
corresponding to the stable fluctuations and described by the coordinates q n (t). The "zero 
mode" (translational degree of freedom) is clamped and frozen since we are only interested 
in describing the dynamics of the unstable mode. The interaction of the stable degrees 
of freedom with the collective coordinate describing the departure from the critical radius 
will introduce viscosity effects and slow down the growth of a supercritical bubble. This 
will result in a smaller growth rate and our aim is to precisely compute this viscosity effect 
on the growth rate for small departures from the critical radius. By including the finite 
temperature static counterterms in the Lagrangian and requesting that these be cancelled 
consistently in the perturbative expansion, we can isolate the static renormalization to the 
unstable frequency from the dynamical viscosity effects associated with the energy transfer 
to the stable modes. 

This program begins by expanding the field in terms of the unstable and stable coordi- 
nates (after clamping the translational mode) as 

<t>(x,t)=<t> b (x,R c )+q- 1 (t)U- 1 (x)+Y,Q P (t)ty>(x), (2.28) 

v 

where the summation index p runs over all stable, bound and scattering, states other than 
the translational and unstable modes. Throughout the index p refers to both scattering and 
bound states. 

Using the above expansion, eq. Q2.28|) , of the field (f>(x,t), the Lagrangian can be shown 



to have the following form: 

L[q-i,Qp] = L [q-i}+ L [q p ]+ L^q^^p] (2.29) 

where 

Lo[q-i] = \ {<f-i(t) + ^ 2 q-i(t) + 5tfq\(t) + hq^(t)} (2.30) 

LM = \ll {q P (t)q- P (t) - " 2 P q P (t)q-p(t)} (2.31) 
/ p 

Li[q-i,q p ] = -V\ - q-i^) \ ^*52B P Qp(t) + X ^V?p 

p,p' p p,p' 

-qt^t) jv^B-i + X^Apq p (t)j- XA- 1 qt 1 (t)+ h.o.t + L ct [q_ 1 ,q p ], (2.32) 

where we have used the equations of motion satisfied by the bubble configurations, and 
equations ( |2.17[ ) and ( 2.18|) . The terms 5Q 2 q 2 i 1 ; hq_i are the quadratic and linear terms 



in q_i from the counterterm Lagrangian. The terms denoted by h.o.t. in ( |2.32| ) correspond 
to higher order interactions that are not important to the order that we are studying. The 
term L ct [q_i,q p ] arise from the non-linear (cubic and higher) terms in the coordinates from 
the counterterm Lagrangian, they do not need to be specified since they will be requested 
to cancel tadpole and static terms in the equations of motion for the collective coordinate 
g_i and do not contribute to lowest order (O(X)). The matrix elements are given by 
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1 f+oo 

B pp , = -j= dx' U_ x {x')U p {x')U v \x')V''\^ ', R c )), (2.33) 
2v A J-co 

1 f + OO 

B p = —= dx 1 U 2 _ 1 (x , )U p (x')V , "((f) b ^\ Rc)), (2.34) 
2v A J-co 

1 f+oo 

B_ 1 = _jf_^ dx'^^F'"^^',^)), (2.35) 

/+oo 
dx 1 U\{x')U v {x')U v i{x'), (2.36) 
-oo 

/+oo 
dx' U\{x')U v {x'), (2.37) 
-oo 

/■+0O 

A-i = / ete' Wl^a/), (2.38) 

«/ — oo 

so that the A-dependence is explicit in each term of the effective action. 

We thus see that the above Lagrangian describes a quantum mechanical degree of freedom 
q-i interacting with a bath of infinitely many degrees of freedom q p , as well as with self- 
interaction. The self interaction of the coordinate g_i will be neglected because we are only 
interested in the viscosity effects arising from interaction with the stable degrees of freedom. 
The non-linear terms in q_i are neglegted because we will extract the growth rate for small 
departures from the critical bubble. 

The bath of mesons will introduce viscosity effects on evolution of the coordinate associ- 
ated with departures from the critical radius and will result in a correction to the growth rate. 
These effects can be studied in a consistent manner by integrating out the meson degrees of 
freedom in a consistent perturbative expansion, thus obtaining the non-equilibrium effective 
action for q_\. The variational derivative of this effective action leads to the equation of 
motion that includes the effects of the meson bath. This non-equilibrium effective action is 
obtained up to one loop order in appendix B. The effective equation of motion is a Langevin 
equation with a non-Markovian viscosity kernel and a Gaussian noise term (to one loop 
order) whose correlations are related to the viscosity kernel via the Fluctuation-Dissipation 
relation. The details are provided in appendix B. 

Alternatively we obtain the equation of motion for the expectation value of the coordinate 
g_x directly by means of linear response. 

We have already accounted for the static renormalization of the unstable frequency Q by 
introducing the finite temperature counterterms and using the finite temperature effective 
potential in the equations of motion. However, we anticipate that there will arise further 
corrections to the growth rate Q from velocity dependent terms in the effective equations of 
motion. These are truly dynamical viscosity effects which cannot be captured by a static 
calculation. The viscosity terms that are a function of the time derivative of the unstable 
coordinate will be revealed by obtaining the effective non-equilibrium equation of motion for 
this coordinate by integrating out the stable modes, which act as a bath for the quantum 
mechanical degree of freedom We will restrict our study to small departures from the 
critical radius and obtain the linearized equation of motion for the expectation value of 
g_i. This is consistent with the definition of the growth rate as the logarithmic derivative 
of the radius for small departures from the critical value. As is clear from ( [2.30D , in the 
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absence of interactions the Hamiltonian for the unstable coordinate g_i is that of an inverted 
harmonic oscillator (of negative frequency squared) and small departures from the unstable 
equilibrium value q_i = will grow exponentially with the growth rate fl. 



F. The initial value problem: 

We study the time evolution of the expectation value of the unstable coordinate Q = (q~i) 
by proposing an initial state described by a density matrix which is a tensor product of the 
density matrix for the unstable coordinate and a density matrix that describes all the stable 
modes in thermal equilibrium at a given temperature. Therefore 

p(0) = p_ 1 (0)®p s (0) (2.39) 

Where p_i describes a pure state in which the expectation value of the unstable coordinate 
is Qo, and p s (0) is a thermal density matrix for the stable modes. As mentioned above the 
translational mode is clamped. The time evolution of the initially prepared density matrix 
is given by Liouville's equation, whose solution is 

p(t) = U(t)p(0)U' 1 (t) (2.40) 

with U(t) being the unitary time evolution operator. As described in detail in reference 
f36| the time evolution can be cast in terms of a time dependent Hamiltonian in which the 
interaction is switched on at the initial time t = 0. Alternatively the initial value problem 
can be cast in terms of linear response to an external source adiabatically turned-on from 
t = — oo that determines the initial preparation 57]. Both approaches are equivalent and 



the reader is referred to for more details on the initial value problem. 

The equation of motion in real time is obtained by using the generating functional of 
non-equilibrium Green's functions which requires a path integral along a contour in complex 



time and the following effective Lagrangian [pq .p 



L eff = L[q± lt q;]-L[qZ 1 ,q;] (2.41) 

where the labels ± refer to the forward (+) and backward (— ) branches along the complex 
time contour [^,^]. The equation of motion for the expectation value Q(t) = (g_i(t)) is 
obtained by performing the shift 

q%(t) = Q(t) + q ± (t) ; (q ± (t)) = (2.42) 

Imposing the condition {^(t)) = to all orders in perturbation theory leads to the retarded 
equation of motion for Q(t). For a detailed presentation of this method in many other 
situations see [ p5[| . 

The important ingredients in this program are the real time Green's functions for the 
stable coordinates which are assumed to be in thermal equilibrium. These are the following 
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%{t)y{t')) = -i5 p ,_ p ,G--(t,t') 

q P (t)q P >(t')) = i5 p ^ pl G-+(t,t') 



-iS. 
-i5. 



p-p 



p-p 



g>(t,t')e(t-t') + g<(t,t')e(t'-t) 
g>(t,t')Q(t'-t) + g<(t,t')e(t-t') 



-i5 p - p ig<(t,t') 
-i5 p ^ p ,g>{t,t') 



$p-p' gp{t ; ,t), 



(2.43) 



where 



2uo r , 



2u] r , 



1 + n p ) exp {— iu p (t — t')} + n p exp {iu p (t — t')} 

'1 + n p ) exp {iujpit — t')} + n p exp {— iuj p {t — t')} 
1 



n,. 



exp {(3u p } 



(2.44) 



We will carry our derivation of the equation of motion in the linear theory where we 
will neglect nonlinear terms, 0(Q 2 ) and higher, and furthermore we will neglect the self- 
interaction of the unstable coordinate (cubic and quartic terms). As mentioned before, 
the linearization of the equation of motion is consistent with the focus on the dissipative 
corrections to the growth rate, which is defined for small departures from the critical radius. 

In this case, the non-equilibrium effective action S e ff becomes 



S, 



eff 



^ dt |L [g + (t)]+L [g p + (t)] + q + (t) 
2\ eW^A^q+Wq+WQit) + h 



-Q(t) + (Q 2 + 5Q 2 )Q(t) - V\ 0(t) £2V<?p + (t)q}(t) 



p,p' 



p,p 



VX9(t)Q(t)J2B pp/ q+(t)q+(t) - [(+) 



p,p 



(2.45) 



where we have written only the terms that are relevant to the lowest order calculation 
O(X) which is the focus of the present discussion and included the time dependence of the 
Hamiltonian to set up the initial value problem [Q. 

Imposing the condition (<?^i(t)) = 0, up to O(X), we obtain the following equation to 
one-loop order 



•oo 
oo 



Q(f ) - n 2 Q(t') - AQ(f) - H - 2i\ B pipA 



p-iPi 



dt(q + (t)q + (t 

' P1,P2,P3,P4 

dt" (<<(0<(0><<(0<(0> - (^(0^(O><ffi(0s£:(O> W) 



X 



with 



A = 5Q 2 - 2\J2A,p> < #(0#(0 > 
p,p' 

n = h- Va^V < £(0<£(0 > 



(2.46) 
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Since the correlation {q + (t)q + (f)) does not vanish for all values of t and t', this implies that 
the quantity between the square brackets must vanish |39 . Furthermore, we now choose the 



counterterm in such a way to ensure that 7i = 0. We thus obtain the equation of motion 
for the expectation value of the unstable coordinate for t > 

Q(t) - Q 2 Q(t) - AQ(t) - A f dt' E(t - t')Q{t') = ; Q(t = 0) = 0; Q(t = 0) = Q 

Jo 

(2.47) 

where the self energy E(£ — £') is given by 
\B I 2 

- O = J2 {(! + n p + v) sin tK + «v)(* - *')] - K - v) sin tK - ^vX* - 01} > ( 2 - 4 ^ 



/-'■/-' 



and we have used the fact that B* , = B- p - p >. The two different terms in the self en- 
ergy, proportional to the sum and difference of frequencies respectively have a simple but 
important interpretation. The first term proportional to the sum of frequencies corresponds 
to the process in which the coordinate g_i "decays" into two mesons with probability 
(1 + n p )(l + n p >) minus the "recombination" process with probability n p n p >. The second 
term proportional to the difference of frequencies originates in Landau damping and cor- 
responds to the scattering of the unstable coordinate with mesons in the medium, with 
probability (1 + n p >)n p minus the reverse process with probability (1 + n p )n p / [flDf] . 

The counterterm SQ 2 will be chosen to cancel the tadpole contribution to A (see eqn. 
( p.46| )) as well as the static contribution from the self-energy, since this contribution is a 
static renormalization of the unstable frequency associated with a stationary bubble. 



G. Viscosity corrections to the growth rate 

In order to obtain the influence of thermal and quantum fluctuations on the growth rate 
of the bubble, we must solve the equation of motion (|2.47 ). This is achieved via the Laplace 
transform. Introducing the Laplace transforms of Q(t), £(£) as Q(s), E(s) respectively with 
s the Laplace variable. The Laplace transform of ( |2.47|) with the specified initial conditions 
is given by 

Q( s ) = s Jk (2.49) 

V ; s 2 - (Q 2 + A) - AE(s) V ; 

where E(s) is given by 

m = £ i^L 2 j (1 + n + ry) 2 /; + 7 r- 2 - K - V ) "j~ Up ' ) (2.50) 

^ u p u p/ { s 2 + (u p + UJ p/ ) 2 " S 2 + {0Jp - Up,) 2 J 

We now isolate the static contribution by substracting 

E(0) = lim s ^t(s) (2.51) 
from the Laplace transform of the self-energy. Using the identity 
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we now fix the counterterm 5Q 2 such that A + AE(0) = 0. Introducing 

f {s) = y\^ll 1 + n p + n p' f! u _ n yp I 1 ) s l 

^ U p UJ pl \ Up + Up* S 2 + (Up + Up,) 2 ' P \Up - Up,) S 2 + (U p - Up') 2 

(2.53) 

the Laplace transform of the equation of motion becomes 

« w - ^_«£W (2 ' 54) 



In order to avoid cluttering of notation we will not write explicitly the V in eqn. ( [2.53| ) 



in what follows, but it must always be understood that the term 1/ (lu p — ay) actually refers 
to its principal part. This principal part prescription arises from the subtraction in the limit 
of vanishing s which is the equivalent of the frequency in the real time domain, since the 
Laplace transform requires that s — > iu + e with uj the frequency . 

At this stage it becomes clear that the procedure of absorbing the local (zero frequency 
limit) contributions to the self energy in the static renormalization of the growth rate Q 
provides the correct description of the dynamics. Viscosity and dissipative effects only arise 
from the time dependence of the expectation value of the coordinate and the static medium 
renormalization had already been absorbed into the definition of the growth rate as the limit 
of zero frequency. Viscosity and dissipation arise from the transfer of energy of the bubble 
wall to other excitations (mesons) and the growth slows down because of these processes. An 
important distinction arises in this case as compared to the familiar situation in field theory 
in which complex poles determine the renormalization to the mass from the real part of the 
self-energy on shell and the decay width from the imaginary part of the self-energy on shell. 
In this case, however, because the frequency associated with the unstable mode is imaginary, 
the real part of the self-energy will renormalize the growth rate whereas the imaginary part 
(if any) at the position of the (purely imaginary pole) will provide oscillatory contributions 
to the bubble dynamics. Therefore, we emphasize that viscosity effects that will diminish 
the growth rate are determined by the real part of the self-energy at the position of the 
pole. This is a striking difference from the usual case in which damping and viscosity are 
associated with the imaginary part of the self-energy on shell. 

The real time dynamics of the bubble growth, Q(t), is given by the inverse Laplace 
transform 

Q(t) = — f e st Q{s) ds , (2.55) 



2iri Jc 



where Q(s) is given by eq. ( |2.54| ) and C refers to the Bromwich contour running along the 



imaginary axis to the right of all the singularities of Q(s) in the complex s-plane. 

The analytic structure of Q(s) consists of cuts along the imaginary axis in the s-plane 
and poles. The two different processes of decay into meson pairs (and recombination) and 
Landau damping discussed above yield two different cut structures. 



18 



The first term in eq. ( |2.53| ) gives a two-meson cut and the contribution from Landau 



damping determines a cut structure that includes the origin in the s-plane. The structure 
of these cuts depend on the matrix elements of the interaction as well as the full spectrum 
of excitations and will be investigated in detail for particular cases in the next sections. 

The growth rate of the bubble with the quantum fluctuation effects included is given by 
the pole s p of Q(s), eq. (|2.54|) , which satisfies the following relation 



s 2 p - VI 2 + XT(s p ) = 0. (2.56) 

To first order in A, the pole of Q(s) which corresponds to a growing bubble is given to 
lowest order (0(A)) by 



s 



v 



±(n + xn {1) ) (2.57) 



where the first order quantum and thermal fluctuation correction Af^ 1 ) to the bubble growth 
is given by 



Afi« 




--£ 

2 , 'jj ri uj r ,' 
p,p' p p 



LU p — tOpi \ Vt 2 + {ujp — UJpi \ 



2 



^2.58) 



We note that both terms inside bracket in ( |2.58j ) are positive since n p < n p i for uo p > ui p > , 
we recall again that the term 1/(uj p — uy) = V(l/ (uj p — uj p i)) from the discussion following 
eqn. ( |2.53|) . Therefore we conclude that the correction to the growth rate is negative i.e. the 



dissipative effects of the coupling to mesons results in a smaller growth rate of supercritical 
bubbles. 

Neither the frequencies of oscillations around the stationary bubble nor the matrix ele- 
ments are quantum mechanical in origin. However the one loop contribution O(X) to the 
self-energy is of quantum origin. 

If we restore h in the expression above this results in A — > Xh ; T —>■ T/h. Obviously 
neither the eigenvalues of the fluctuation operator nor the matrix elements depend on h but 
the one loop contribution to the self-energy includes the h from the coupling (one loop) as 
well as from the temperature factors. 



H. Classical Limit 

In the next sections we will see that the correction to the growth rate is dominated by 
low lying excitations and for high temperatures these will be such that u p « T. In this 
case we can invoke the classical limit which is best understood by restoring h as mentioned 
above 

1 T 

n k = — sa - — >> 1 (2.59) 
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In this form 



n k -n kl ^ — [ ) = ~T ( 2 - 60 ) 

h \u k uj k i ) ft u k uj k i 

and in the high temperature limit we further approximate 

1 + n k + n k , « 1 + Tl ±±^L „ ^^±^ (2 . 61) 

These approximations lead to a simplified expression for the lowest order correction to the 
growth rate 

^ = _^y lVI 2 J ^ + ^ \ (2 62) 

2 ^ c^, \ O 2 + + c^ P 2 ^ 2 + - u; p ,) 2 J ' 

where the product AT is independent of ft and displays clearly the classical limit. 

The classical limit will be justified in each particular case in the next sections. The 
expression ( p. 62] ) is the final form of the one loop 0(A) correction to the growth rate arising 
from dynamical viscosity since all of the static contributions had been absorbed by the 
counterterms. This is as far as we can pursue in a general manner without addressing the 
details of the spectrum of fluctuations around the bubble solution. In the next two sections 
we study the details of the model determined by the Lagrangian ( |2.1| ) with the potential 
721) for the cases of 1 + 1 and 3 + 1 dimensions. 



III. THE 1+1 DIMENSIONAL CASE 



As mentioned in the introduction, the 1 + 1 dimensional case is relevant in statistical 
and condensed matter physics. Quantum field theory models based on the Lagrangian ( |2.1|) 
with potentials with a stable and metastable state are proposed to describe the low energy 
phenomenology of quasi-one dimensional charge density wave systems [O], and therefore 



their relevance in these physical situations warrants the study of this case. 

For V{4>) given by eq. (|2.2|) , the solution to the static classical equation of motion 
for one spatial dimension can be found exactly [ |I3| . |H| 
given by |M| 



The critical bubble is found to be 



%(x, s ) = + 



m 



2V2A 



tanh 



+ s 



— tanh 



s 



2 

m' 



(3.1) 



where £ is the width of the bubble wall and Sq is given in terms of the critical radius R c by 



s = 



R c 

T 



-cosh 1 



1 



(3.2) 



with e given by eq. Q2.5|) . The above solution corresponds to a kink-antikink pair centered 
at x = 0, and separated by a distance R c , and is displayed in fig. (g). This is the one 
dimensional bubble that starts at the false vacuum d>_ almost reaches the true vacuum d>+ 
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and returns to 0_ and it is similar to the polaron solution found in quasi-one- dimensional 



polymers [T3| . 

The total energy of the bubble, given by eq. fl2.11| ), as a function of its dimensionless 
radius s can be calculated from the field 4>b(x, s ) by replacing s = R c /(, —> s = R/£ in eq. 

The gradient term (d(j)b/dx) 2 contributes to the surface energy while V (</>&) has two con- 
tributions; surface contribution and volume contribution. Substituting V (</>&) in eq. ( J2.ll ) 
and evaluating the integral, one finds that the total energy of the bubble is given by 



E„ 



E sur (s) + E vol (s) 



with 



—Am 3 



4 m s 



+ 



3m 3 (-1 +e 8 



8e 4s s) (2 + 77) 



3 A (-l + e 4s ) A ' 4(-l + e 4 -) z A 

4e 4s m 3 (l-e 4s + 2 (l + e 4s ) s) 



+ 



-l + e 4s ) 3 A 



E vo i(s) = 
where r\ is defined as 



rj m 



(3.3) 



(3.4) 
(3.5) 



n 



2. 



(3.6) 



The volume contribution grows linearly with the radius of the bubble while the surface 
contribution saturates at about the critical radius of the bubble, Sq, and attains the following 
asymptotic value 



lim E sur (s) 

s^oo 



m 3 (2 + 9 rj) 



12 A 



In the thin wall limit (77 — > 0), the surface energy is simply twice the kink mass [p5|l , as the 
kink-antikink pair becames widely separated and the exponential interaction between kink 
and antikink becomes negligible. 

The total energy of the bubble is depicted in fig.(|3]). It attains its maximum, E c , at the 
critical radius so where 



m 



24A 



4 + 12 77 + 3 rj 2 - 3 s V (§ + 6 ij + rj 2 ) 



with 



so 



-cosh 1 
2 





"4 + r/" 




77 



(3.7) 



(3.8) 



In the thin wall limit, the above expressions simplify and we find 

fe-D 2 



V 



+ 0((e 
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and 



E r 



m 
6A 



+ 0((e 



(3.9) 



To study the fluctuations around the bubble solution we need the complete set of eigen- 
functions {U n } satisfying 



d 2 



dx 2 



+ V"( 



UJx) 



io n U n {x) 



(3.10) 



with 



V"{ 



3 m 2 



m 



-sech z 



so 



3 m 2 



-sech z 



s 



Although the spectrum of eigenfunctions and eigenvalues is known exactly in the case of 
one kink or antikink |35 |, for the case of the kink-antikink pair there are not known results 
that we are aware of. Solving for the eigenfunctions {U n } in this case is a difficult task but 
in the thin wall limit, the above potential consists of two identical and widely separated 
wells centered at x = ±i? c . The spectrum of each potential well is known in the literature 
|35| . It consists of a zero frequency mode localized in the well, an excited bound state with 



frequency 3m 2 /4 that is also localized in the well and a continuum of scattering states. 

We can use approximate methods such as the linear combination of atomic orbitals 
approximation (LCAO) (available in elementary textbooks) to provide a reliable estimate 
for low lying bound states of the above potential from the spectrum of the single potential 
well. In this method, the ground state of the potential V"((pb) is the symmetric linear 
combination of the ground states of the single well located at x — ±R C while the first 
excited state of V"((pb) is the anti-symmetric linear combination. 

In the thin wall limit, these two states are given by 



U-i(x,s ) 



and 
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2^/E~ c ds 



K x,s ) 



with 
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-n 2 



-24mV 4s ° = -6r]m 2 (3.11) 



Uo(x, s ) 



r Er dx 



(x,s ) 



with 



0. 



(3.12) 



These are the unstable, U-i, and the zero, U , modes which we discussed earlier, see eqs. 
( j2.19| ) and ( |2.21|) . Obviously U Q (x) is associated with translations since it is the spatial 
derivative of the bubble solution and must correspond to a vanishing eigenvalue by trans- 
lational invariance. Since it is antisymmetric there has to be a nodeless eigenfunction of 
smaller frequency. The symmetric combination is U-\(x) and since the two combinations 
will be split off in energy by a tunneling amplitude that is exponentially small in the dis- 
tance between the kink and the antikink, the unstable frequency must be negative and 
exponentially small in this separation as is clearly displayed in eqn. ( |3.11| ). 

Besides the unstable and the zero modes, there are two bound states that have energies 
~ (3m 2 /4) ± AE(sq) with energy difference which is again exponentially small in s an d 
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correspond to the symmetric and antisymmetric linear combinations of the bound states of 
the kink in 4 theory ||3f| localized at x = ±i? c . Since V'^b) does not exactly reach m 2 near 



the center of the bubble, the potential in the Schroedinger equation could allow for shallow 
bound states near the scattering continuum with binding energies that are exponentially 
small in the variable sq. These bound states if present are extremely difficult to obtain. 

Finally there is the scattering continuum region of the spectrum characterized by func- 
tions Uk, whose eigenfunctions are asymptotically phase shifted plane waves with eigenvalues 
u\ = k 2 + m 2 

In order to compute the matrix elements that enter in the expression for the correction 
to the growth rate ( |2.62| ) we need either the exact form of the eigenfunctions or an excellent 
approximation to them. Whereas we are confident of our analysis regarding the low lying 
bound states, we lack a full understanding of shallow bound and continuum states. Such 
an understanding requires a detailed study of the spectrum which certainly lies beyond the 
scope of this article. Although we do not have a complete understanding of the spectrum 
of eigenfunctions and therefore we do not have even a good approximation to the matrix 
elements, we can, however, provide some physically reasonable assumptions complemented 
with dimensional arguments to provide an estimate for the corrections in this case. 

We begin by noting that in the one kink case, the potential that enters in the Schroedinger 
equation for the fluctuation is reflectionless |H| and the scattering states for the one kink 
(or antikink) case have a transmission amplitude which is a pure phase. 

Thus we expect that in the thin wall limit when the kink-antikink pair separation is 
much larger than the width of an isolated kink the wave functions of the scattering states will 
acquire a phase shift that is at least twice as large as that in the one kink case and will have a 
rather smooth dependence on the kink-antikink separation. Because the potentials for each 
kink are reflectionless we assume that in the thin wall limit the reflection coefficients of 



each potential well are very small and therefore there are no substantial interference effects in 
the region between the kink and the antikink. Under these assumptions the matrix elements 



B P)P i will be similar to those calculated in reference |41| for a single kink case and fall off 
very fast at large momenta justifying the classical limit [41|. Hence, under these suitable 
assumptions the matrix elements are smooth functions of Sq = R c /£,. The contribution from 
the two-meson cut, i.e. the first term in the bracket in ( |2.62|) is proportional to Q because 



the f2 in the denominator can be neglected in comparison with the frequencies for the meson 
states uj p ~ 0(m). Since in 1 + 1 dimensions the coupling A has dimensions of (mass) 2 the 
correction to the growth rate arising from the two meson cut is of the form 

AT 

Sn 2mes « MF[s ] (3.13) 

with F[sq] a dimensionless slowly varying function of So = R c /£, which is rather difficult to 
calculate and can only be obtained from a detailed knowledge of the eigenfunctions. 

The contribution from the Landau damping cut is more complicated to extract. The 
second term in the bracket in ( |2.62| ) has the form of a Lorentzian and since Q is exponentially 



small it is a function that is strongly peaked at u v = u p > and the sum (integral) over p, p' is 
dominated by a region of width Q near u p = u p >. Assuming that the matrix elements are 
smooth functions of momentum, in this one dimensional case the integral over a small region 
ojp — Up' ~ Q can be done by taking a narrow Lorentzian and integrating over the relative 
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momentum within this region [42], leading to a contribution of the form — (XT/m 2 )C[so\ 
where C[sq] is a smooth function of its argument that can only be calculated from a detailed 
knowledge of the scattering wavefunctions. Hence 



5Q 



LD 



XT 



(3.14) 



and the total shift in the frequency is given by 5Q = 5Vt2 mes + SQld- 

Thus we conclude that although we do not have a complete knowledge of the eigenfunc- 
tions and therefore cannot provide a complete calculation of the correction to the growth 
rate, suitable assumptions based on the properties of the spectrum of the one kink case 
combined with dimensional arguments suggest that to lowest order in the coupling and in 
the thin wall approximation the growth rate of a slightly supercritical bubble is given by 



n 



4^ _Mc 



8 1 f J 



(3.15) 



with i* 1 ^], C[^] slowly varying dimensionless functions of their argument. Obviously a full 
calculation of F[-y],C[-y] requires a detailed understanding of the spectrum, in particular 



the scattering states. However it is clear from (|3 .15 ) that the validity of a perturbative 
expansion places a severe constrain on the coupling constant and the value of the temper- 
ature, in particular the second contribution in ( 3 . 1 5|) arising from Landau damping gives 
the leading correction in the thin wall approximation and signals a potential breakdown of 
perturbation theory in this limit. A more detailed understanding of this possiblity requires 
a better knowledge of the scattering matrix elements, this is an extremely difficult problem 
that depends on the details of the potential and lies outside the scope of this article. 



IV. THE 3+1 DIMENSIONAL CASE 
A. General aspects: 

We now study the 3 + 1 dimensional case which is more relevant from the point of view 
of particle physics, however before focusing on a particular form of the potential and bubble 
profile, we can study fundamental model independent properties of the 3 + 1 dimensional 
case that will determine very robust predictions for the corrections to the bubble growth. 

The static bubble configuration 0j (r, R c ) is radially symmetric and satisfies the static 
equation (|2.10| ). 

To study quantum fluctuations around the critical bubble configuration, we need to find 
the spectrum of the fluctuation operator M. which in 3 spatial dimensions is given by 



M = -V z + 



2 , d 2 V((f>) 



(4.1) 

<f>b(r,Rc) 



Since the critical bubble solution is radially symmetric, we write the eigenf unctions 
Mnim( r , 9, V 9 ) °f the differential operator M. as a product of spherical harmonics Y [m (9,ip) 
and radial functions ip n i{r) that satisfy 
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dr 2 



2d_ 1(1 + 1) d 2 V{ct> b ) 
r dr r 2 dd> 2 



Mr) = J^niir) 



(4.2) 



Because of the translational invariance of the Lagrangian ( |2.1|) , there is a 3-fold degen- 
erate zero mode given by the eigenfunction oc V0&(a?, R c ) oc Yi^^d^/dr which correspond 
to translations of the bubble in 3-dimensions with no energy cost. These are the Goldstone 
modes of the spontaneously broken translational invariance. This can be easily seen by 
taking the derivative of eq. ( p.lOj ) with respect to r which results in the following equation 



cP 



dr 2 



2d_ 2_ d 2 V{<t> b ) 
r dr r 2 dd> 2 



dr 







(4.3) 



In addition to the Goldstone mode, there are other low-lying excitations corresponding 
to different values of I ^ 1. The eigenfunctions of these excitations and their eigenvalues 
can be obtained by writing the / = 1 term in the above equation as 



1(1 + 1) (/-!)(/ + 2) 



and we re-write equation ( |4.3| ) in the following form 



dr 2 



2d_ 1(1 + 1) <9V(0 b ) 
r dr r 2 dd> 2 



dr 



{i -m+ 2) 

R 2 



a<Pb + S y,, acpb 
dr dr 



where 



5V(r) = (/-!)(/ + 2) 



1 



(4.4) 



(4.5) 



(4.6) 



Since the function dfo/dr is strongly localized at r = R c in the thin wall limit, the second 
term on the right hand side is a small localized perturbation. Therefore the unperturbed 
lowest lying eigenvalues are given by 



01 



(*-!)(* + 2) 
Rl 



(4.7) 



In the thin wall limit, £/R c << 1 we find that the lowest order correction (in SV) to 
these eigenvalues is of order 0(^ 2 /R 2 ). For details see fT4] , fl"3"| and appendix This analysis 
reveals that there is a band of low lying modes with eigenfunctions 



Uoi m (r,9,ip) = JN iY lm (6,Lp)- 



(r,Rc) 



dr 



(4.8) 



with N Q i the normalization constants, corresponding to the 2/ + 1-fold degenerate eigenvalues 
given by ( |4.7|) . Using equations (|2.16|J2.26| , p.21| ) and (|4.7|) with I = (corresponding to the 
unstable mode) we find the normalization to be given by 
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Rio 



(4.9) 



in terms of the critical radius and the surface tension. 

We summarize several noteworthy features of these low-lying solutions: 
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These excitations become Goldstone modes in the limit as R c — > oo, i.e. in the limit 
in which the radius of curvature of the bubble goes to infinity. This statement will 
be understood in detail below in connection with the case of flat interfaces in 3 + 1 
dimensions. 

The eigenfunction with the lowest eigenvalue, corresponds to I = 0, i.e. a spherically 
symmetric solution with a negative eigenvalue given by |14|,|3],|43| 

ul, = -n 2 = - A (4.io) 

This is the unstable mode which corresponds to a spherically symmetric expansion 
or contraction of the bubble and therefore corresponds to the unstable functional 
direction. The coordinate associated with this mode is the displacement from the 
critical radius. 

The three Goldstone modes ZYo,i,o ; ZYo,i,±i are the translation modes. 

The higher energy modes with I > 2 are excitations on the surface of the bubble, or 
surface waves with energies given by 

^= (/ ~ 1 ^ + 2) =y(^-i)(^ + 2); l>2 (4.11) 



These low lying modes will play a dominant role and we will refer to them collectively as 
Uoim with eigenvalues given by ( (4.7| ) and whose normalized eigenf unctions ( |4.8p are simple 
functions of the bubble configuration which for the potential (|2.2|) in the thin wall approxi- 
mation are given by eqn. ( f4.27|) below. 



That the modes with / > 2 can be identified as wiggles of the bubble surface, or sur- 
face waves, can be seen from the following expansion [|43] , fi4|l where as discussed before the 
translational modes is not included because it is "clamped" 



(r,t) = <f>b(r- R c ) +q-i(t)Uooo(r,0,(p) + J2 a im{t)Yi m (6,(p) — 

l>2;m ar 



+ ■■■ 

Rc 

Mr - R(9, <p, t)) + g_!(t)W oo(r, 9,<p) + --- (4.12) 



where 



R(e,(p,t)=R c - Yl ai m (t)Y lm (e,(p) 

l>2;m 

This is an important identification that we emphasize: this band of low-lying modes 
describes fluctuations of the surface of the bubble. We will argue below that these excita- 
tions dominate the infrared behavior of the viscosity correction and will provide the largest 
contribution to the viscosity coefficient. 

The low lying spectrum described above is fairly general and only depends on the exis- 
tence of a thin wall bubble. Depending on the form of potential V((f>), there might be other 
bound states. 
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An analysis similar to that leading to the band of low-lying excitations reveals that for 
the potential ( j2.2|) there is another band of rotational bound state excitations that starts 



4 



3m 2 /4. In the thin wall limit the radial wave function for the lowest state in this 



band is given by the bound state of energy 3m 2 / '4 of the 4 theory in 1 + 1 dimensions |35 



which is also localized at the wall. The eigenfunctions and eigenvalues for this rotational 
band of excitations are given by (see appendix A) 

U llm (r, 6, <p) = JNuY lm (e, p)Mr) ! <4 « ^ + + 0(4~) (4.13) 



4 m "r, 



c 



where ipi(r) is given by the bound state of the theory with potential (|2.2|) in one spatial 
dimension |[35|| , the eigenfunctions in this band are given in eq. (|4.28|) below. 

Finally there is a continuum of scattering eigenstates with eigenvalues io\ = k 2 + m 2 . 

As will be discussed later, the contribution to the growth rate from the rotational band 
( |4.13j ) and of the scattering states is sub leading in the thin wall limit. 

The maximum value of angular momentum available for the low-lying part of the spec- 



trum (|4.7|) is limited by the edge of the continuum spectrum or the presence of higher bound 
states, hence l 2 max lR 2 c < rn 2 or l 2 max < {mR c ) 2 = (i? c /£) 2 . Therefore in the thin wall limit 
the maximum value of the angular momentum l max » 1. 



B. Planar interfaces, surface waves and (quasi) Goldstone bosons 

The low lying spectrum of eigenfunctions given by (|4.g| ) with eigenvalues ( |4.7|) has a 
simple physical origin, which can be understood by noticing that in the limit R c — > oo the 
discrete spectrum becomes a continuum. 

In the limit when the radius of the critical bubble is very large R c — > oo, the interface 
between the two phases + and 0_ becomes planar and the two phases become degenerate. 

Let us consider a static planar interface configuration corresponding to a domain wall 
along the z-axis in three spatial dimensions. Such a configuration satisfies the following 



equation ||35,43 



-^# + ^M = (4.14) 

dz 2 ocp 

where (f) w (z) satisfies the boundary condition 4> w {z) — > <f>± as z — > =Foo and z is the coordinate 
perpendicular to the planar interface. 

The quantum fluctuations around the classical static wall solution 4> w (z) are given by 
the spectrum of the differential operator 

M = -V* + *y (4.15) 

Since the domain wall only depends on the coordinate z, the differential operator for 
the fluctuations is separable in terms of eigenfunctions ^^(z, x±) = e^ J -' x±, ipn(z), where x± 
denotes the transverse coordinates to z, namely x and y, and q± the transverse momentum. 
The functions ip n {z) are so lution of the following eigenvalue problem 
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d 2 



dz 2 



d 2 V(<j) w (z)) 



(4.16) 



Taking the derivative of eq. (|4.14|) with respect to z, i.e. 





d 2 d 2 V(^ 



dz 2 



dz 



(4.17) 



and comparing to eq. (|4.15|) we see that d(f> w /dz is the zero mode which corresponds to the 



translational invariance [35|. 
Therefore the eigenf unctions 



Vto = exp (iq ± ■ x ± ) 



dz 



(4.18) 



have eigenvalues q\. These are Goldstone modes associated with translational invariance 
and represent excitations of the surface of the planar interface <fi w (z) since 



<t>{r) = <p w ( z ~ z o) + J2 S± ex P (iq±-x±) 



dz 



~ <p w (z - z±(x±)) 
where z is the position of the planar interface and 



z ± (x ± 



20 



q± 



a q± exp (tq±.x_ 



This is clearly similar to the case of the spherical bubble eq. ( 4.12| ) and describes the same 
physics, i.e. fluctuations of the surface. In the case of a planar interface these surface waves 
are also called capillary waves, and describe the hydrodynamic modes of long- wavelength 
fluctuations of interfaces in systems with two degenerate phases separated by an interface 
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In the case of degenerate phases, such as for example a liquid-gas or an Ising system, 
these surface waves are Goldstone modes associated with the breakdown of translational 
invariance by the presence of the interface. 

For a spherical bubble in the thin wall limit these surface waves acquire a gap given 
by the inverse radius (proportional to the Gaussian curvature of the surface). Therefore in 
analogy with the case of interfaces for degenerate separated phases, we identify these surface 
fluctuations as quasi Goldstone modes. Since, as argued before the maximum frequency of 
the surface waves is < m they are classical in the high temperature limit T » m. Hence 
the surface waves are identified as classical hydrodynamic fluctuations of the bubble shape 
and quasi- Goldstone modes in the thin wall limit. 

We want to emphasize that these low energy excitations are a robust feature of the thin 
wall approximation and are model independent. Having studied in detail the general aspects 
of fluctuations around a thin wall bubble, we now focus on the specific details of the theory 
with a potential given by ( |2.2j ) so as to be able to compute the matrix elements and provide 
a quantitative analysis of the viscosity effects. 
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C. 4> A theory: specifics 



The critical bubble solution that satisfies eq. ( |2.1Up with the potential ( |2.2| ) is not in 
general an elementary function, but in the thin wall limit the critical radius of the bubble 
mR c ~ R c /£ >> 1 and the function dfo/dr is localized near R c which makes the "friction" 
term, r~ 1 d<j)b/dr oc l/(mR c ) ~ £/R c << 1- in this limit, the critical bubble solution is 
found to be 



k(r, R c 
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(4.19) 



It corresponds to a field configuration that starts around r = at the true vacuum + , 
given by eq. (|2.6| ), and goes to the false vacuum, 0_, as r — > oo with a surface width 
£ = m/2 and a critical radius R c . Having specified the critical bubble solution <p b , we now 
go back and determine explicitly the general expressions which we discussed in the previous 
sections. 

The total energy of the bubble as a function of its radius R, given by eq. (|2.11| ), can 
be calculated from the field 4>b(r, R c ) by replacing R c in eq. flOp by R. The gradient 
term (V(/>&) 2 contributes to the surface energy while V(4>b) has two contributions: a surface 
contribution K (</>&) an d volume contribution K>(0&) given by 
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(4.20) 
(4.21) 



where rj is defined by eq. (|3.6| ). Substituting the above expressions in eq. (|2.11| ) and 
evaluating the integral, we find that the total energy of the bubble is given by 



= — V (((>+) R 3 + 



i + _ ^1 

2 ~ 2 i? 



V-6 37T 2 ?/ \ ^ 
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(4.22) 



where in the thin wall limit 



V(<f>+) 



rj m 



+ 0((e 



The total energy attains its maximum at 

1 



Rr 



2 + 9 rj + J 4 + 36 rj + 45 rf 



(4.23) 



Using the fact that rj is small in the thin wall limit, we find that the critical radius R c is 
given by 



Rc 



3^ 



l + 0((e-l) 



(4.24) 



and the total critical energy is given by 
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(4.25) 



which is equivalent to eqs. (|2.13|) and (|2.14 ) with the surface tension given by 

,3 



a 



m 
12A 



(4.26) 



The low-lying fluctuation modes £Y z m (#, y?, r) eqn. ( f4.8|) are given by 

r-R c 



Uoim(r, 8, (p) = 6 m sech 2 
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and the next rotational band of bound states is given by (for details see appendix A) 



Uu m (r,8,ip) ~ J Nu sech 
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(4.28) 



D. Corrections to the bubble growth rate 

From eq. (|2.62p , the corrections from quantum and thermal fluctuations to the bubble 
growth in the present case has the following form in the classical limit 
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Since V" is spherically symmetric, the angular integral leads to 
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(4.30) 



(4.31) 



(4.32) 



M r )%>i{r) (4.33) 



We now note that the bound state energies only depend on /, and that the energy of the 
scattering states is independent of I. Therefore there is no Landau damping contribution 
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from the bound states as consequence of the principal part prescription which subtracts 
the contribution from uo p = uy in the Landau damping term as discussed in detail below 
equations ( |2.53| , 2.58|) . Therefore the correction to the growth rate can be written as 5Q = 



Sflf, + 8Q S with the correction from the bound states given by 

AO V y (2/ + l)|ffnx,*| 2 n , . 

b ~ 2 1 ^ <n W + K* + ^) 2 [ } 

and that from the scattering states given by 

The bound state correction has three contributions; the quasi- Goldstone modes contri- 
bution corresponding to surface waves with = f2 2 (Z — 1)(Z + l)/2; Z > 2, the contribution 
from the higher energy rotational band of bound states near u; 2 > 3m 2 /4 and the mixed 
contribution from the quasi- Goldstone and the higher energy bound states modes. The de- 
nominators in 5Qb are of order Q 6 for the quasi-Goldstone modes contribution as compared 
to m 6 for the mixed and the higher energy bound states contributions. Since in the thin wall 
approximation m/Q oc R c /£ » 1, the largest contribution arises from the quasi-Goldstone 
surface modes with the lowest energy denominators. It can be easily seen that the matrix 
elements cannot compensate for the difference in powers of R c and that in fact for the higher 
energy bound states these matrix elements are smaller than those for the surface modes be- 
cause the wave functions Uu m actually vanish at the position of the bubble wall (see eqn. 

a ). 



The contribution from the scattering states is also seen to be much smaller than that 
from the surface waves. The frequencies for the continuum states uj^ > m » Q and the two 
meson cut give a contribution of order Q (since in the denominators the Q can be neglected 
as compared to m). For the Landau damping cut, the argument is similar to that of the 
case of one space dimension. 

For Q << m this contribution has a Lorentzian shape of width m Q, and the integral 
over the momenta can be performed in the narrow width approximation. The subtraction 
of the static contribution guarantees that the integral is dominated by the Lorentzian W2 



and the region uo^ w uiy can be integrated by changing to relative variables and now in 
three spatial dimensions the phase space in the region Uk — toy ~ f2 give extra powers of 
Q as compared to the one dimensional case. Furthermore, the matrix elements are smooth 
functions of the radius of the bubble as can be understood simply by a scattering argument 
from a sharply peaked potential in three spatial dimensions. These matrix elements do not 
introduce any singularity in the limit R c /£ » 1, hence the contribution of the scattering 
states is at least proportional to Q and is therefore subleading in the thin wall limit. 

This analysis leads to the conclusion that the largest contribution to the correction to 
the growth rate is given by the quasi-Goldstone modes, i.e. the surface waves, since these 
are the lowest lying excitations and hence provide the smallest energy denominators. 

The matrix element Bo o i for the surface waves can be calculated easily and we find 



vr 5R 2 C 
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leading to the following correction 
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(4.37) 
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where 



21 + 5 



(Z + 1) 2 (/ + 4) 2 (1 + 2(Z + !)(/ + 4)) 



~ 0.039 



(4.38) 



1=0 



The above series converges rapidly and only the first few terms contribute to the sum, 
we have evaluated the sum numerically with R c /£ = 10. 

Hence we summarize one of the main results of this article: the lowest order correction 
to the bubble growth for three dimensional bubbles is dominated by viscosity effects arising 
from the excitation of long-wavelength surface waves and is given by 



Therefore to lowest order in A (one-loop) and to leading order in the thin wall approxi- 
mation we find that the growth rate of slightly supercritical bubbles is given by 



with £ and R c are the width and radius of the critical bubble. This is one of the important 
results of our study. The validity of perturbation theory places a very stringent constraint 
on the quartic coupling constant in the thin wall limit i? c /£ >> 1 and in the classical 
limit T/m ~ T£ >> 1. The validity of the classical limit in this case is warranted: we 
are studying the dynamics of nucleation via thermal activation for temperatures below the 
critical temperature T < T c ps m/y/X but for temperatures much larger than the energy of 
the low lying excitations, the relevant regime for thermal activation is m/y/X > T >> m in 
the weak coupling limit. Furthermore the low- lying excitations with frequencies « m are 
obviously classical. 



The focus of this article is to provide a microscopic calculation of the growth rate of 
slightly supercritical nucleation bubbles. The model under consideration is a <f> A scalar theory 
with an explicitly symmetry breaking term that produces a metastable and a stable ground 
state, we studied the case of nucleation in 1 + 1 dimensions as well as 3 + 1 dimensions. 
The former is relevant in the case of quasi-one dimensional charge density wave systems 
and organic conductors. We begin our analysis by obtaining the critical bubble solution by 
including finite temperature effects in the potential that enters in the classical equation of 
motion, counterterms are added to the Lagrangian to compensate for the finite temperature 
corrections consistently in a perturbative expansion. 




(4.39) 




(4.40) 



V. CONCLUSIONS AND DISCUSSION 
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Our approach to obtaining the growth rate is very different from previous treatments in 
that we begin by expanding the quantum field around the critical bubble in terms of the 
quadratic fluctuations around the critical bubble configuration. These fluctuations describe 
an unstable direction associated with small departures from the critical radius, translational 
zero modes and stable fluctuations. The translational modes are anchored by fixing the cen- 
ter of the bubbles and we treat explicitly the interaction between the coordinate associated 
with the growth (or collapse) of the bubble with those associated with the stable fluctua- 
tions. We obtain the growth rate by obtaining the effective linearized equation of motion for 
the unstable coordinate by integrating out the coordinates associated with the stable fluc- 
tuations. Two different approximations are involved; i) a weak coupling expansion in terms 
of A the quartic self-coupling, and ii) the thin wall approximation in terms of £/R c with £ 
the width of the bubble wall and R c the critical radius. The first approximation allows a 
consistent perturbative expansion of the self-energy of the unstable coordinate, the second 
allows a quantitative calculation of the relevant matrix elements, furthermore our analysis 
reveals that the important fluctuations are classical for temperatures T c > T » m, with 
T c the critical temperature and m the mass of quanta in the metastable phase. 

In the one dimensional case we are able to provide an estimate for the growth rate given 
by eq. ( |3.15| ) where the functions F[R c /£] , C[R c /£\ depend in a detailed manner upon 
the scattering states solutions of the eigenvalue problem for the quadratic fluctuations, and 
clearly will depend on the details of the potential. This estimate points out the potential 
breakdown of perturbation theory in the thin wall limit. 

In the case of three dimensions we are able to extract some robust features that trascend 
the form of the potential and are solely a consequence of the thin wall limit. In particular we 
identify a rotational band of low lying excitations which describe surface waves i.e. ripples on 
the surface of the bubble. We establish the connection of these surface waves to the capillary 
waves of flat interfaces in the case of degenerate but phase separated thermodynamic states 
(such as the Ising or liquid - gas at coexistence), the surface waves are then identified as 
classical long-wavelength hydrodynamic fluctuations. The unstable coordinate couples to 
these hydrodynamic fluctuations and as a result the friction term arising from the self- 
energy is dominated by the coupling to these hydrodynamic modes. Clearly the coupling 
of the unstable coordinate to these hydrodynamic modes depends on the model, and in the 
case under consideration we find in the thin wall limit and to lowest order in perturbation 
theory the following expression for the growth rate 



We also obtain the effective non-Markovian Langevin equation for the coordinate de- 
scribing small departures from the critical radius and establish the generalized fluctuation 
dissipation relation between the viscosity and the noise kernel. The noise is correlated on 
time scales comparable to Q^ 1 precisely as a consequence of the coupling to the hydrody- 
namic modes and cannot be treated simply with white (delta function) correlations. 

Discussion: Although we have studied a specific microscopic model, in the case of 3 + 1 
dimensions we have been able to identify some robust features that trascend the particular 
model. These are the existence and dominance of hydrodynamic fluctuations associated 
with surface waves in the thin wall limit. The coupling of the coordinate associated with 




(5.1) 
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small departures from the critical radius to these low energy fluctuations induces friction 
or viscosity corrections to the growth rate of slightly supercritical bubbles and in the weak 
coupling and thin wall limit these fluctuations give the largest contribution to the friction 
corrections. 

One of our original motivations is to make contact with previous studies of nucleation as 
applied to the quark-hadron phase transition. In particular Csernai and Kapusta || have 
parametrized the coarse grained free energy that describes a quark-hadron first order phase 
transition in terms of a local energy variable that can be identified with our scalar field </>. 
The form of the potential taken by these authors coincides with our potential V(<p) ( f2.2p . 
with coefficients that depend on temperature, just as we have argued in this article. Their 
form of the critical bubble solution and the variational energy as a function of the radius 
of a bubble are very similar to those studied in this article. These authors have established 
that for about 1% supercooling the critical radius for such a model is about R c ~ 12fm the 
width of the wall is about £ ~ 0.7fm and the thin wall approximation must be reliable in 
this regime. 

We can obtain the quartic self-coupling A in ( |2.2| ) from the parameters used in || by 
relating the width of the bubble £ = 2/m and the surface tension a to A via eqn. ( |4.26|) . In 
H the value of the surface tension for the particular quark-hadron model is a = 50MeV/fm 2 
for T w 200MeV, yielding a value A ~ 6 which leads to a very large correction to the 
growth rate. Certainly the large value of the coupling invalidates the perturbative scheme 
and we cannot draw a definite conclusion as to the relevance of our lowest order estimate 
for the quark-hadron transition beyond the statement that the viscosity induced by the 
hydrodynamic fluctuations could result in a very large (negative) correction to the growth 
rate and therefore a rather small nucleation rate. 

For larger supercooling the critical radius becomes smaller and the thin wall approxima- 
tion breaks down, but in this case nucleation and spinodal decomposition will be indistin- 
guishable and homogenous nucleation theory may not be the proper description. 

However despite the limitations of the perturbative expansion and the thin wall approx- 
imation, we have provided a consistent approach to obtain friction or viscosity corrections 
to the growth rate from a microscopic perspective without invoking a phenomenological de- 
scription. The observation that in the thin wall limit the most important corrections arise 
from the coupling to classical hydrodynamic fluctuations could perhaps pave the way to 
a systematic hydrodynamic treatment that would circumvent the weak coupling expansion 
and allow to extract non-perturbative physics. 
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APPENDIX A: CORRECTIONS TO THE QUASI-GOLDSTONE MODES 



In this appendix we calculate the corrections to the eigenvalues of the quasi- Goldstone 
modes that obey the eq. (|4.5|) using first order perturbation theory. First we write 
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First order energy correction Ei are given by 



E, 



(i) 



{U l m \5V\U i m ) 

Rl 



dr 



(A2) 



where we have integrated out the angular degrees of freedom. 

For the 4 potential given by eq. ( |2.2| ), the normalized wavefunctions are given by eq. 
( [4.271 ) an d when it is substituted in eq. ( |A2| ) one finds that 
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Thus to first order in perturbation theory, the eigenvalues are given by 

(/-!)(/ + 2) 



(A3) 



0l>, 



01 



Rl 



l + O 



Rl 



A similar treatment can be used for the next rotational band based on the bound state 
of the 4 theory in one spatial dimension. For this we treat the term 
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r dr 
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(A4) 



as a perturbation. The unperturbed wave function is the positive energy bound state of the 
one dimensional <j) 4 theory [35] and the first order correction is obtained as before and is 
l(l + l)/Rl + 0(-L). 



APPENDIX B: LANGEVIN EQUATION AND FLUCTUATION DISSIPATION 

RELATION 

The semiclassical Langevin equation is obtained by performing the path integrals over 
the bath degrees of freedom, i.e. the stable modes, thus obtaining a non-equilibrium effective 
functional for the unstable mode. This is achieved consistently in perturbation theory |39| , f4T 
and to lowest order we find 
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Z[j+=r = 0] = / VqZ 1 VqZ 1 e^^i^-Lo^) jrfe+^j, ( B l) 
where the Lagrangians Lo[<7-i] are given by eq. ( |2.30| ) and qZi] is the influence 



functional |^5| , ^6| which to one loop order and neglecting the tadpole contributions is given 

by 

F[q±i, Q-i] = expj - i 2 \ ^IVI*/ dt dt' [g+ (t) t>)G+ + (t, t') <?+(*') 



+ qZ 1 (t) G~ (t, t')G-r (t,f) qZ 1 (0 - <£i (*) G+- (t, t')G + r (t,f) qZ x (f) 



gz 1 (^)G- + (^,t')G p - + (^,t / )^l(^ , ) 



(B2) 



where the Green's functions are given by eq. ( |2.43j ). 

Introducing the Wigner coordinates or center of mass and relative coordinates, x(t) and 
r(t) respectively given by 

<t) = \ + Q-M) > r(t) = qMt) - qZ 1 {t). 



the generating functional becomes |46|]- |50 



Z[0] = J VxVr e iS[x ' r] 
with the non-equilibrium effective action given by 

S[x,r] = J dtr{t) -x(t) + Q 2 x(t) + \J dt' (s ret (t - f) x(t') + % -K(t -t')r(t'] 



(B3) 
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where for clarity we have neglected the counterterms as well as the tadpoles arising in the 
computation of the influence functional. The kernels S re ^(t — t') and K(t — t') are given by 



s ret(*"*') = Z(t - t')Q(t - t') 

K{t-t>) = - £ \B P1P2 \ 2 {g>(t,t')g>(t,t') + g<(t,t')g<(t,t')} 
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with S(t - f) given by (|2l^ ). 

At this stage it proves convenient to introduce the identity 



e 2. 



1 / dtdt'r(t) K(t-t')r(t') _ f e -±Jdtdt'£(t)K- 1 (t-t>)t(t>)+iJdt£,(t)r(t) 



with C(t) being an inessential normalization factor, to cast the non-equilibrium effective 
action of the unstable mode in terms of a stochastic noise variable with a definite probability 



distribution |46fl- [50| . Using the above relation, the generating functional becomes 
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Z[0] = J VxVrV£P[£\ expji J dt r(t) -x(t) + Q 2 x{t) + A J* dt' £(t - t') x(tf) + f(t) |, 

(B6) 

where the probability distribution of the stochastic noise, P[£], is given by 

P[t] = Jvt exp {-^ j dtdt'i(t)K~\t - t'W)} ■ (B7) 
In this approximation we find that the noise-noise correlation function is given by 

«£(t)£(0»= K(t-t), (B8) 

which is in general colored, i.e. it is not a delta function S(t — £'). 

The semiclassical Langevin equation is obtained by extremizing the effective action in 
eq. (|B6|) with respect to r(t) 



x{t) - Vt 2 x{t) - A / dt' '£(i - if) x{t') = f(t) 



(B9) 



Taking the average of the above equation with the noise probability distribution P[£] and 
identifying << x(t) »= Q(t) yields the equation of motion for the expectation value of 
the unstable mode, eq. (|2.47|) . 

The relationship between the kernels £ re t(i — t') and Kit — t') constitutes a generalized 
quantum fluctuation dissipation relation. This relation is established by considering the time 
Fourier transforms of the functions Qp(t, t')Gp>(t, tf) and Qp(t,t f )Q5(t,t') which are denoted 
by Qppiiyj) and Qp P i{yj) respectively. These Fourier transforms obey the KMS condition |51 
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Using the above relation we find that E re ^(u;), the Fourier transform in time of the 
retarded self energy S re ^(t — t') is given by 
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leading to the imaginary part 
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On the other hand the kernel that determines the noise-noise correlation function K(t—t') 
has a Fourier transform given by k(u) with 



fc( W ) = (i+ e -^)x;iv%(w) 

~(3u 



coth 



im[E ret (u;)]. 
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The above relation between the Fourier transform of the noise-noise correlation function and 
the imaginary part of the self-energy is the generalized fluctuation-dissipation relation. 
In particular the contribution from the surface waves to the kernel K(t — t') is given by 

Km(T) = ^ | -J®±}L- ( , + cos [2^m^)£ } (B14, 

This function oscillates with constant amplitude on a time scale ~ R c which is the same 
time scale for growth of a supercritical bubble. Therefore the noise term cannot be taken 
to be uncorrelated over the time scales associated with the growth of a bubble, i.e. the 
noise is colored and a Langevin description based on a white noise would miss the long-time 
correlations. The contribution of the scattering states could lead to a short range part of 
the kernel, but the long time behavior of the kernel will be dominated by the low energy 
surface waves. 

The reason that a Markovian Langevin equation with white noise fails to describe the 
dynamics of the unstable coordinate is that there are slow hydrodynamic fluctuations with 
time scales comparable to the growth rate that couple to the unstable coordinate, precisely 
the same type of fluctuations that dominate the viscosity or friction. 
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FIGURES 
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FIG. 4. The total energy of the three dimensional bubble as a function of its dimensionless 
radius s for e = 1.2, A = 0.1 and 0_ = 2.0. 
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